A temperature behavior of the frustrated translational mode of adsorbate and the 
nature of the "adsorbate— substrate" interaction 
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A temperature behavior of the frustrated translational mode (T-mode) of a light particle, coupled 
by different regimes of ohmicity to the surface, is studied within a formalism of the generalized 
diffusion coefficients. The memory effects of the adsorbate motion are considered to be the main 
reason of the T-mode origin. Numerical calculations yield a thermally induced shift and broadening 
of the T-mode, which is found to be linear in temperature for Ohmic and super-Ohmic systems and 
nonlinear for strongly sub-Ohmic ones. We obtain analytical expressions for the T-mode shift and 
width at weak coupling for the systems with integer "ohmicity" indexes n = 0^2 in zero temperature 
and high temperature limits. We provide an explanation of the experimentally observed blue- or 
rerf-shifts of the T-mode on the basis of a comparative analysis of two typical times of the system 
evolution: a time of decay of the "velocity-velocity" autocorrelation function, and a correlation time 
of the thermal bath random forces. A relation of the T-mode to the multiple jumps of the adsorbate 
is discussed, and generalization of conditions of the multiple hopping to the case of quantum surface 
diffusion is performed. 

PACS numbers: 05.60.-k, 62.25.Jk, 63.22.-m, 66.10.cg 
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I. INTRODUCTION 

Among a number of problems dealt with the surface 
diffusion, which occurs either due to thermahy activated 
[l|, Ql or tunneUing mechanisms [21, i3l| , one can single out 
the investigation of the frustrated translational mode (T- 
mode) of the adsorbate [J, Q • The T-mode is known to 
appear when the adparticle moves between two stable 
positions at the surface within one adsorption site (for 
instance, between the atop position and the saddle point 
[3|) due to interaction with phonons or electronic sub- 
system of the substrate. This motion of the adsorbate 
is being detected during inelastic helium atom scatter- 
ing (IHAS) [y, 01 as an additional peak in the dynamic 
structure factor [a,H,@|. The vibrational nature of the T- 
mode allows to refer this excitation to the external (low 
frequency) modes of the adsorbate [Ifl, [ll|, which are 
more or less temperature dependent and dispersionless 
with homogeneous broadening at low coverages [5J. 

In the last years the T-mode became a subject of 
the intensive studies of both experimentalists and the- 
orists. For instance, in the recent scanning tunnelling 
spectroscopy experiments ^21 the excitation of the T- 
mode in polyatomic molecules CO/Pd(110) has been ex- 
plained in terms of the energy transfer from the C-0 
stretch mode to the frustrated translational one due to 
the anharmonic coupling between normal modes. A cer- 
tain intrigue is introduced by the fact that the T-mode is 
observed even in such "classical" systems as CO/Cu(001) 
and CO/Pt(lll) T3|, though the underlying mechanism 
of the T-mode formation is purely quantum mechanical: 
the transitions between ground and excited states within 
a well. This stimulated a creation of the models based 
on classical [13J or quantum jl4j] Langevin equations, the 
transition state theory (TST) [i^|, where anharmonic- 
ity of the lattice potential is taken into account, leading 



to the temperature induced shift and broadening of the 
T-mode fS]. 

At the same time, a theoretical approach for the T- 
mode description has been proposed [5| from purely 
quantum mechanical concepts. This method allows 
to extrapolate the results obtained at T ^ to the 
zero temperature limit, yielding information about the 
"adsorbate-substrate" interaction strength and vibra- 
tional frequency from the IHAS data. Much effort has 
been put Q into description of the T-mode tempera- 
ture dependence. It is known that unlike diffusion the 
vibrational damping is not an activated process (5|, [l6| . 
Nevertheless, the temperature induced T-peak shifts to- 
wards lower frequencies, and almost linear broadening is 
observed [13 for the Na/Cu(001) system in a region be- 
tween K and 200 K, while for the CO/Cu(001) and 
CO/Pt(lll) systems the T-mode is 6Zue-shifted with 
temperature. It is believed that at low temperatures 
the shape of the T-mode is dominated by a frictional 
damping, characterized by the nonadiabatic coupling to 
the substrate excitations p]. At high temperatures the 
anharmonicity of the static lattice potential plays a dom- 
inant role 4]. Summarizing, it could be stated that 
up to now there is no unique viewpoint which factors 
(phonon and/or e-h excitations, anharmonic contribu- 
tions, or memory effects) are predominant in the T-mode 
development at different temperatures, couplings and in- 
teraction regimes. 

Last but not least, one can relate the existence of the 
T-mode to multiple jumps of the adsorbate. The onset of 
multi-hops was associated with an inelastic peak of the 
dynamic structure factor at frequency uJosc of the adpar- 
ticle oscillations at the well bottom [i|,|l3l- If the T-mode 
is well resolved on the dynamic structure factor and not 
overlapped by the central quasi-elastic peak, it could be 
said that thermalization of the adsorbate velocity has 



not taken place yet, and multiple jumps contribute sig- 
nificantly to the diffusion coefficient. 

In the present paper we study the conditions of the 
T-mode formation, its temperature dependence, and its 
connection with the nature of the "adsorbate-substrate" 
interaction. We explore a two-level model of the light 
particle coupled non-adiabatically to the surface phonons 
|18l42(j |. The other kinds of "adsorbate-substrate" inter- 
action, like electronic friction or nonlinear phonon cou- 
pling, will be also considered. The underbarrier hopping 
to the nearest adsorption sites is postulated to be a "driv- 
ing force" of the surface motion of the adsorbate. 

We consider non-Markovian effects of the adsorbate 
motion, dealt with retarded relaxation of the lattice ex- 
citations, to be the main reason of the T-mode origin. 
Since this one-particle excitation is very sensitive to the 
form of the kinetic kernels, rigorous numerical evalua- 
tion of the T-mode temperature behavior is carried out 
without any approximation (e.g., like Wigner-Weisskopff 
one [3, lU). On the other hand, both zero tempera- 
ture and semiclassical limits for the kinetic kernels are 
considered in our paper that allows us to obtain the ana- 
lytical expressions for the T-mode characteristics at low 
"adsorbate-substrate" coupling. 

We perform our analysis in terms of the generalized 
diffusion coefficients (which are related to the "velocity- 
velocity" autocorrelation functions) rather than dynamic 
structure factor. The latest approach is often com- 
plicated due to the overlap of quasi-elastic and inelas- 
tic peaks of the dynamic structure factor [ij, while the 
method of generalized diffusion coefficients allows one to 
study the T-mode in more detail, being focused on the 
adsorbate motion at the well bottom. We propose a phys- 
ical interpretation of the observed temperature induced 
shift of the T-mode on the basis of a comparative analysis 
of two typical timescales of the system dynamics: a time 
of decay r„ of the generalized diffusion coefficients and a 
correlation time Tcor of the thermal bath random forces. 
The obtained results reveal a close connection between 
the direction of the T-mode shift and the nature of the 
"adsorbate-substrate" interaction. 

Our paper is organized in the following way. In Sec. II 
we present the generalized diffusion coefficient of a light 
particle derived by the method of quantum kinetic equa- 
tions [i^ [231 on the basis of two-level dissipative model 
of adparticle coupled to substrate phonons. In Sec. Ill 
we consider different correlation functions of the thermal 
bath random forces, which depend on the nature of the 
"adsorbate-substrate" interaction. In the next Section 
we attribute the T-mode of the adsorbate to the coher- 
ent part of the generalized diffusion coefficient and calcu- 
late its temperature dependence numerically at different 
values of ohmicity index n. A comparison of the ob- 
tained results with the experimental data and the results 
of similar theoretical approaches is provided in this Sec- 
tion. In Sec. V we calculate the T-mode frequency and 
width analytically in zero and high temperature limits 
for strongly sub-Ohmic (ohmicity index n = 0), Ohmic, 



and the super-Ohmic systems with n = 2. In Sec. VI, 
according to Refs. [l|, [l3l, we relate the T-mode to the 
onset of multiple jumps of the adsorbate and general- 
ize the conditions of the multiple hopping to the case of 
quantum surface diffusion. In the last Section we discuss 
briefly the obtained results and draw final conclusions. 



II. GENERALIZED QUANTUM SURFACE 
DIFFUSION COEFFICIENTS 



We consider a light particle, which performs under- 
barrier hopping to the nearest adsorption site, oscillates 
between ground and excited states within the potential 
well, and interacts with the lattice vibrations. Hamilto- 
nian of the system is chosen as follows [la, [l^ 



H — Ha + Hint + Hb, 



(2.1) 



where the adsorbate is described by the two-band con- 
stituent 



i?A=^(-ioalo«s'0+iialias'i)+^ — (n^i - nsQl2.2) 

{ss') 



Y 



Here s denotes the site of the lattice; and 1 are the 
ground and excited states within a given well, and {ss') 
denotes a sum over the nearest-neighbor sites. The quan- 
tum states within a well are referred to as "vibrational" 
ones with the vibrational frequency 51. a^j (osi) creates 
(destroys) a particle on the site s in the vibrational state 
i; Usi = OsiOsi is the number operator for this state, and 
Us = Uso + risi. Hereafter we will deal with a single 
adparticle only, hence the adparticle statistics becomes 
irrelevant, io and ti are the nearest-neighbor tunnelling 
amplitudes in the ground and the first excited states, re- 
spectively, and we expect that ii ^ io- 

The coupling to phonons is considered to be local 
within each well. Phonons may couple both to the ad- 
sorbate density operators and to the vibrations within a 
quantum well. The interaction Hamiltonian is ,18.] 
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Isqibq + bl) + {al^asi + al^aso) 



E^^^(^^ + ^J) 



(2.3) 



where 6j (bq) creates (destroys) a phonon with a nor- 
mal mode frequency coq. The strengths jsq iXsq) de- 
scribe coupling of phonons to the density (oscillation) 
modes of the adsorbate. The bandwidths to, ti and vi- 
brational frequency J7 can be evaluated in the framework 
of the eigenvector-eigenvalue problem for a periodic po- 
tential, felt by an adsorbate due to the static lattice. 
The coupling strengths are expressed via the mean values 
r — {s,i\Vfnt\s,j) of the lattice distortion potential Vj^nt 
over the localized Wannier states \s,j) times the phase 



factor depending on the site number s and wave-vector q 
[18| . Likewise, we suppose T to be the same for different 
quantum states {i,j} = {0, 1} and use the dimensionless 
couphng parameter 



G = 



r2 



MujI 



(2.4) 



to characterize the "adsorbate-substrate" interaction. 
Here uj,nax stands for the Debye frequency, and M de- 
notes the mass of the substrate atom. 
The last term in Eq. (|2?T|) 

(2.5) 



Hb ^^ fuJqblbq 



corresponds to the phonon bath; longitudinal acoustic 
phonons only are taken into account in this model. 

Since the tunnelling amplitudes are always much 
smaller than the coupling strengths, the interaction part 
of Hamiltonian p.3p cannot be considered as a perturba- 
tion and has to be taken into account exactly. The most 
refined method consists in performing unitary transfor- 
mations in order to exclude the linear over the interac- 
tion terms from the system Hamiltonian. In the unitary 
transformed Hamiltonian the tunnelling processes can be 
considered as those with emmision/adsorption of the vir- 
tual phonons, and the same is true for the vibrational 
processes. 

It is also useful to pass to the hybrid set of states for 
each site: 

a^L=—{asQ±asi), (2.6) 

and similarly for the creation operators. The designation 
L or R means that a single adparticle is now localized on 
the left or right side of the given well. We will refer to 
the vibrational transitions with i ^ j, {i, j} ~ {L, R}, as 
the end-changing processes, supplying them afterwards 
by the subscript (c), and transitions with i — j will be 
termed as the end-preserving ones with the corresponding 
subscript (p). 

Using the method of the reduced density matrix [21| 
it is possible [la, [l^, l23 to obtain the chain of quantum 
kinetic equations for diagonal fs,s{t) — X]i=L fl('^L'^s«)s 

and off-diagonal fs,s'{t) = ^i^i, Ri^l'i'^siYs one-particle 
non-equilibrium distribution functions, where the aver- 
aging is taken with the statistical operator ps{t) of the 
adsorbate. These integro-differential equations are lin- 
ear in /s^s(t), fs.s'{t) but are non-local in time; hence it 
is useful to perform the Laplace transformation f(z) = 
L exp{—zt)f{t)dt. Solving the equations for the off- 
diagonal distribution functions and inserting the ob- 
tained results into the equation for the diagonal ones, 
one can obtain |19l . |20[ an expression for the generalized 
(frequency dependent) diffusion coefhcient: 



D{z) = D,oh{z) 
4 



f A. 
2*2 






Z + linteA'A -K ^intrai'A + ^thi'A 



Winter \^ X 



..7) 



The first term Dcohiz) in (|2.7p describes a coherent con- 
tribution to t he g eneralized diffusion coefficient. It can 
be interpreted [la,[l3 '^^ terms of a simple model of band- 
type motion limited by scattering from the lattice (with 
interatomic spacing a) at temperatures large relative to 
the bandwidth tinter , which is narrowed due to polaronic 
effect. The kinetic kernel 



linter[z) = A^Ll{z) + 2^lr{z) + 2^rl{z) 



(2.8) 



corresponds to the dissipative intersite motion of the ad- 
sorbate and describes processes, when the adparticle per- 
forms random site-to-site hopping (with or without the 
change of its quantum state) owing to the interaction 
with the bath. The kinetic kernel jmtraiz) in Eq. (|2.7p 
describes a dissipative intrasite dynamics, when the ad- 
sorbate during its scattering from the lattice gets enough 
energy from the bath to be excited from the ground state 
to the upper level within the same adsorption site (the 
opposite process of particle de-excitation with a phonon 
emission is also taken into consideration). 

The second term Din{z) in Eq. (|2.7p is an inco- 
herent contribution to the generalized diffusion coeffi- 
cient. This is the result expected from the random 
walk model for diffusion with site-to-site hopping rate 
linter{z), describing processes of the surface phonon cre- 
ation/annihilation, when the particle performs a transi- 
tion from one Wannier state to another. 

The rates jintra{z), "/inter (z) can be obtained from the 
Laplace transformation of the kinetic kernels 



7a;(r) =: a;„iaa;A2Re{exp[-((^j;(0)-(^3;(T))]-exp[ 



(2.9) 



iLir)-- 



where 



fxir)-^ 



U}maxtlRe{exp[-{ipLLiO) + ^LLiT))]-exp[~ipLLiO)]}, 



Jxiuj) 



coth 



2k^T 



(2.10) 



cos(a;T) — i sin(wT) (?.ll) 



and one-particle parameters Aa; = {^1,17} are related 
to the corresponding end-ch ang ing /end-preserving pro- 
cesses (see also Table in Ref. [101 )■ The exponential form 
of the rates (|2.9p - (|2.10p is a result of averaging of the 
lattice time correlation functions [l6J, [l9[ over the bath 
variables that yields the temperature dependent factor 
in Eq. (|2.1ip . Hereafter we use dimensionless frequencies 
in the units of ujmax and temperatures in the units of 

fli^max/kB- 

The "adsorbate-substrate" interaction enters the ki- 
netic rates via spectral weight functions [la, Il8l - [20| 

>/M = E^'/('^-^?)' (2.12) 

JhRii^)^^ [{isq - ls'q) + {Xsq + Xs'q)f5{uJ - W,), 



(2.13) 



jRL{i^)=^[insq ~ls'q)~{Xsq+Xs'q)] 5{uj-0Jq), 



JLLii^)^"^ Klsq - ls'q) + {Xsq " Xs'q)fS{uj ~ UJq), 



(2.14) 



Jrr{^) =X! [(T«9 ~ ls'q)~{Xsq " Xs'q)f5{u} - Wq). 



The function (|2.12p describes the intrasite dynamics; the 
functions ()2.13p are related to the intersite end-changing 
processes, while (|2.14p are dealt with the intersite end- 
preserving processes. The spectral weight function (I2.12p 
can be considered as site-independent if the system has 
a translational symmetry, whereas the functions (|2.13p - 
(I2.14P depend only on the distance \s — s' \ between the 
nearest neighbor sites s and s' [la, [20| • 

The system dynamics is governed by the low frequency 
behavior of the spectral weight functions (|2.12l) - (|2.14l) . In 
Refs. [ig, Il8l - [20l | they were chosen to be scaled as 



Jc{lo) w e(a; - a;o)9(a;,„ax - a;)r/ca; 



D~2 



(2.15) 



for the end-changing, and 

Jp(a;) sa 6(w - ujQ)<d{uj.rnax -~ ^)'nc^^ ■, (2-16) 

for the end-preserving processes with r^c = lOG, rjp — 
12. 5G, and 8(2:) denoting the Heaviside step function. 

In Eqs. (|2.15p - (l2.16|) the power index D equals the 
system dimensionality. Thus, for a two-dimensional lat- 
tice the end-changing spectral weight functions are sub- 
Ohmic, while the end-preserving ones are super-Ohmic 
[221 ■ Such a behavior is a consequence of the equality 
of the coupling strengths T for all the quantum states. 
Moreover, the lattice is assumed to possess a nonzero 
lowest frequency wq, which is introduced to take into ac- 
count the finite size of the system. It not only removes 
the divergencies [22| . when the sub-Ohmic spectral func- 
tions have the power index n = 0, but also allows one 
to describe the adsorbate-induced surface reconstruction 
[18| . when the particles become self-trapped due to the 
overlap of lattice distortions. 

We would like to note that a sharp cut-off of the spec- 
tral weight functions at the Debye frequency uimax is not 
a unique one, and plain exponential [22| or algebraic [23| 
cut-offs are also used. However, the obtained results are 
found to be rather insensitive to the cut-off form, being 
governed above all by the low frequency behavior of J(w). 

In the next Section we will consider other values of 
the power index n in the expressions for spectral weight 
functions and discuss their relation to the processes of 
different physical nature. Here we would like to empha- 
size that there is a one-to-one correspondence between 
the low frequency behavior of J{uj) and the long-time re- 
laxation of the kinetic kernels (|2.9p . It was pointed out 
in Ref. [20] that at large times we pass from a Gaus- 
sian form for 7a; (t) (for the sub-Ohmic processes with 
n = 0) through exponential relaxation of kinetic kernels 
(for the Ohmic processes) to the power law behavior (in 
the super-Ohmic case with n = 2). 



To conclude this Section, let us note that at a reason- 
able assumption of an infinitesimal value of the tunnel- 
ling amplitude in comparison with vibrational frequency, 
ti ^ Ml, the expression (|2.7p for the generalized diffusion 
coefficient can be simplified and presented as follows 



D{z) 



2+2 r 



a^t 



2h^ 



exp(-y)p(0)) 

z + n^%{z) 



+ 2(7,(z)+7p(z) 



47) 



where ^c{z) and 'Jpiz) mean the Laplace transforms of 
the end-changing (evaluated with spectral weight func- 
tions (J2.15I) ) and end-preserving (evaluated with spec- 
tral weight functions (|2.16p ) kernels (12. 9p times factor 
l/wmaxA^. A generic form for the coherent contribution 
to D{z) has a structure similar to Eq. (|2.17p even if the 
other kinds of interaction (like electronic friction and/or 
anharmonic coupling to phonon subsystem) are taken 
into account. In the very general case, a denominator 
in the expression for Dcoh{z) will consist of several terms 
of different physical origin, which describe the vibrational 
transitions of the adsorbate within the well. In Sec. IV we 
attribute the T-mode to the coherent part of the general- 
ized diffusion coefficient, and explore Eq. (I2.17P to study 
in detail its features at different types of the "adsorbate- 
substrate" interaction. 



III. CORRELATION FUNCTIONS OF THE 
THERMAL BATH RANDOM FORCES 

We can look at the problem of quantum surface diffu- 
sion from another viewpoint. Making use of the solution 
of the Heisenberg equations of motion for the external 
degrees of freedom [2J] one can derive the so-called gen- 
eralized quantum Langevin equation [^, [2^ 



Mq{t) + M / 7(t - T)q{T)dT 
to 



dViq,t) 
dq 



Frit) (3.1) 



for the particle with the mass M and a generalized co- 
ordinate q{t) moving within the potential V{q,t), where 
the damping kernel j{t) is expressed via spectral weight 
function in the following way 



7(t) 



J(c.) 



cos{u!t)duj. 



(3.2) 



The mean value of the operator- valued noise Fr{t) is zero 

{Fr{t))bath^O, (3.3) 

whereas the time correlation function of the thermal bath 
random forces can be expressed as follows 

1 

2' 



SpFit) 



M 

■K 



-,{Frit)Fr{0) + FriO)Frit))i,ath 



-hu coth 



2kBT 



cos{ujt)duj. (3.4) 



In the classical limit it converts to the non-Markovian 
Einstein relation 



S^Mt) = MkBT-/{t), 



(3.5) 



where the average in Eqs. p.3|) - (|3.4p is taken with re- 
spect to a bath density matrix, which contains shifted 
oscillators 25 1. 

The relation p.4p implies that modelling of quantum 
dissipation in the system is possible in terms of macro- 
scopic quantities such as the friction kernel ^(t) and tem- 
perature T. A class of such systems is by no means 
restricted to the "adparticle - linear Debye phonons" 
coupling. One can give examples of several types of 
the "adsorbate-substrate" coupling, which are closely re- 
lated to the low frequency behavior of the spectral weight 
functions. They range from the systems with a flicker 
noise [2^ [23] with J{uj) ^ w° (like the end-changing 
spectral functions (I2.15p ) through the Ohmic coupling 
with J{uj) ~ uj^, when one takes into account the elec- 
tronic friction and/or anharmonic terms in "adsorbate- 
substrate" interaction [2l|,[23|, to the super-Ohmic sys- 
tems with J{uj) ^ up' , when coupling to surface Debye 
phonons is considered [18, 2Q] (see also the end-preserving 
spectral functions (|2.16p ). 

At any ohmicity it is useful to introduce the correlation 
time Tcor which defines the decay rate of correlation func- 
tions p.4p - p.5p of the thermal bath random forces. In 
other words, at times about Taor the excess energy of the 
lattice relaxes. For instance, in the systems with flicker 
noise, the corresponding time correlation functions p.5p 
have the residual correlations Tcor ^ w^^^ in damping 
kernels. 



7(i) -- -T]cCi{u!Qt), 



cost 



Ci(z)= I dt— 



(3.6) 



The sub-Ohmic case p.6p with the power index n = 
implies an introduction of the lower limit cjq of frequency, 
< Wo ^ uJmax, which is related to the inverse time of 
the experiment duration [27| , or dealt with the finiteness 
of the system size [18 2^ (see also Eqs. ((2T5l) - ((2T6| and 
subsequent explanation). 

The Ohmic systems have a white-noise-like behavior 
7(t) ^ rjc5{t), and there are ultra-short range time corre- 
lations in the super-Ohmic case. The super-Ohmic cou- 
pling with the power indexes n > 2 is beyond our consid- 
eration. The cases of non-integer n, which are typical for 
fractal systems, will be considered solely for the reason 
of proper convergence of the numerical calculations. In 
particular, in the next Section we consider the ohmicity 
indexes n = e, and n = 2 — e with e <C 1. An alternative 
way is to use the lower bound cut-off frequency wq, as it 
is done in Sec. V, where the analytical expressions for the 
T-mode features in the quasi-classical limit are obtained. 

There is also another typical time scale of the system 
dynamics - the relaxation time r„ of the velocity auto- 
correlation functions. This time will be introduced in the 



next Section, and along with T^or will be used for expla- 
nation of the direction of the temperature induced shift 
of the T-modc on the basis of a comparative analysis 



between 



and 



IV. T-MODE AS A COHERENT PART OF THE 
GENERALIZED DIFFUSION COEFFICIENT 

As has been already said, the T-mode of the adsorbate 
is being detected in IHAS experiments as an additional 
(inelastic) peak in the dynamic structure factor. Thus, 
to predict theoretically the T-mode features one has to 
evaluate the dynamic structure factor 



^(k,^) = 7^ / exp-''^*/(k,i)di 
27r J 



1 

271 



— I e-*'^*(e-^'"^(°)e''^-^(*))dt. 



(4.1) 



which is nothing but the Fourier transform of the in- 
termediate scattering function /(k, i). In Ref. Q the 
model of a damped quantum harmonic oscillator bilin- 
early coupled to a bath of lattice oscillators has been 
proposed. Non-Markovian equations of motion for the 
creation/anihilation operators of the adsorbate have been 
derived, which have a form of the generalized Langevin 
equation (|3.ip . The dynamic structure factor of the ad- 
sorbate has been found with taking into consideration 
the memory effects. Temperature dependences of the fre- 
quency of the inelastic peak of 5'(k, w) and its full width 
on half maximum (FWHM) have been studied in detail. 

Though an investigation of the T-mode features using 
the dynamic structure factor is closely related to scatter- 
ing experiments, this approach can sometimes be com- 
plicated due to the overlap of quasi-elastic and inelas- 
tic peaks. It happens, for instance fl^, in the quasi- 
continuous regime of adsorbate motion, at small fric- 
tion and high values of /c, when the quasi-elastic peak of 
^(k, w), dealt with diffusion of the adsorbate, becomes 
low and broad and cannot be well separated from the in- 
elastic peak at the vibrational frequency 51. The inelastic 
peak also becomes indistinguishable at strong coupling 
due to the impossibility of completing coherent oscilla- 
tions in the wells. 

In such a case, investigation of the T-mode 
features using the Fourier transform C„t,(aj) = 

l/27rRe J_ e''^*Ct,^(i)(ii of the velocity autocorrela- 
tion function C^„ (i) looks more promising, being focused 
on the adsorbate motion at the bottom of the potential 
well. For instance, C^t,(cj) reveals the adsorbate motion 
in the potential wells also in the overdamped regime at 
high barriers, when, in spite of large friction, the particle 
can perform small parts of oscillations with unthermal- 
ized velocity T]. On the other hand, Cyv{uj) can be easily 



obtained from the dynamic structure factor by a simple 
relation [171] 



t'n^(cj) 



2,. S(k,uj) 
u) lim 



fe->0 



k2 



(4.2) 



It has been pointed out in our previous papers [19|, |20| 
that the generalized diffusion coefficients are directly re- 
lated to the velocity autocorrelation function Cyy{t), de- 
termined at the adsorption site s. We will follow this rule, 
having dealt a time-dependent diffusion coefficient D{t) 
with Cvv (t) . Such a definition of the generalized diffusion 
coefficient differs from that used in Ref. j2S], where the 
relation D{t) — J^ Cyy{u)du has been adopted. However, 
it is much more useful since it allows one to associate the 
time of decay of D{t) with the relaxation time t^ of the 
velocity autocorrelation function. On the other hand, 
the zero frequency limit D^xp — Trlimi^^o C't)i;(w) = 
L D(t)dt determines the experimentally measured dif- 
fusion coefficients Dg^p, whose temperature dependence 
has been studied quite profoundly [2|, |3| . 

In Ref. [20| it has been shown that there is a close 
relation (at least, within the model proposed) between 
the adparticle dynamics at intermediate times t <^ Ty 
and temperature dependence of the diffusion coefficients. 
Namely, as the coupling strength increases, the adparti- 
cle motion (initially oscillatory) becomes more and more 
smooth, indicating that the temperature behavior of the 
diffusion coefficients Dexp{T) should change from weakly 
dependent on T to quite a sensitive function of temper- 
ature. The oscillatory behavior of the generalized diffii- 
sion coefficient has been related to the T-mode onset. In 
the frequency representation, oscillations of Cyy (t) yield 
the side wings of Re[Ci,^(aj)], localized in the vicinity 
of the vibrational frequency fl. The T-mode frequency 
was found to increase with temperature (see Fig. 3 in 
Ref. 20]). The corresponding side peak of Re[Ci,t,(a;)] 
should be Wue-shifted too. 

Thus, there are two peaks of different physical origin in 
the Fourier transform Ke[Cyy{Lu)] of the velocity autocor- 
relation function, presented in Fig. [TJ The central one, 
localized at w = 0, is related to the processes of adsorbate 
site-to-site random hopping. This peak is determined by 
the incoherent contribution 



A„(i) = {ati/h)^Re [-fc{t) + -fp{t)] 



(4.3) 



to the generalized diffusion coefficient. In a weak coup- 
ling limit, at hight temperatures, and at large times the 
end-changing kinetic kernels decay as 

7c(t) ^ exp [-r/el \nujo\ {ksTT^ + ir)] , (4.4) 

while the end-preserving ones decay as 

7p(r) ^ Qiu:,' - t) (l/r^".^-^ - J,^'"^'') . (4.5) 

It is seen from Eqs. (|4.4l) - (|4.5p that the relaxation law of 
the kinetic kernels changes from the Gaussian behavior 
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FIG. 1: Frequency dependence of the Fourier transform 
of velocity autocorrelation function at Debye temperature 
Td = hiUmax/kB, ohmiclty index n — 0, and parameters 

tl/hWmax = 10~^, U>o/u>max = 10"", 77c = 10"^, Q./u}max = 

l/\/2. 



for the sub-Ohmic case with n = to the truncated long 
tails for the super-Ohmic case with n = 2. It is also easy 
to show that the kinetic kernels in the Ohmic regime 
behave as 



7c (r) -- exp(-7r?7cTr). 



(4.6) 



In a weak coupling limit and at the vibrational fre- 
quency large enough there is also an additional maximum 
of Re[C't,^(a;)]. This peak is determined by the adsor- 
bate scattering from the substrate atoms [16] and corre- 
sponds to recrossing processes dealt with the adparticle 
nonmonotonic dissipative motion [l9|, [20| . Hence, a side 
peak arises due to the coherent contribution to the gener- 
alized surface diffusion coefficient (|2.17|) . The expression 
for frequency dependence of the T-mode, associated with 
the above mentioned maximum of Re[Ct,t,(aj)] can be eas- 
ily obtained from Eq. (|2.17p . putting z = —iuj + 0+ 



C,T(^) = Re D,oh{uj) 



2j.2 



a^t 



2R{lo) 



m^ [w-172/(^)]2+f^4^(^)2' 



(4.7) 



/(a;) = / 7c(t) sin(a;i)(it, -R(w) = / "/dt) cos{ujt)dt. 





It is seen from Eq. (14.71) that the T-mode localization in 
the vicinity of the vibrational frequency fl is determined 
by I{lu), whereas the T-mode width depends on R{uj). 

Let us study the time dependence of the coherent con- 
tribution to the generalized diffusion coefficient in more 
detail. Performing the inverse Laplace transformation. 



one can obtain the following expression: 



Dcoh{t) =Rc 



e-\-ioo 



— -— -5— hm / dzexp[zt) — , ^„. , , 



atj 



Re 



^exp(zii) — 



1 



n^^',(z. 



The summation in (|4.8p in accordance with the residue 
theorem runs over all poles Zi of the integrand, which 
obey the condition Re[zi] < 0. In a general case, the 
summation is extended to the infinite number of poles, 
and the major contribution comes from terms with maxi- 
mal values of Re[2:i] and weight factors [1 + 7^(2^)]"^. 

The expression for Dcoh{t) can be even more compli- 
cated if one deals with higher order poles. Indeed, when 
the temperature increases, the 2-nd order poles (with 
rather small weight factors) appear in the integrand. 
Even though the high order contributions to Dcoh{t) have 
a form similar to anharmonic corrections to the velocity 
autocorrelation function [j] , their physical nature is def- 
initely different. Basically, they are of non-Markovian 
origin, while a skewness of the T-mode peak in Ref. [4| 
arises due to the anharmonicity of the local lattice po- 
tential. 

Just as relaxation times of the incoherent contributions 
Din(t) to the generalized diffusion coefficients depend on 
the ohmicity of the system (see Eqs. (|4.3p - (|4.6p ). the 
damping rates | Re Zi | in (|4.8p also change over the large 
range of values depending on nature of the "adsorbate- 
substrate" interaction. An additional study shows that 
at fixed parameters of the model there is a strong in- 
equality, Max{Re[zi(n = 0)]} < Max{Re[2;i(n = 1)]} < 
Max{ Re[zi{n = 2)]}, Re[zi(n)] < 0, ordering the damp- 
ing rates of the coherent term Dcoh{t) subject to the 
ohmicity index n. 

Turning back to Fig. [1] let us point out that the T- 
mode width, which defines the inverse life-time of this 
one-particle excitation, can increase for two different rea- 
sons. At fixed temperature the T-peak is broadened when 
the "adsorbate-substrate" coupling increases, and the co- 
herent motion of the adparticle diminishes. Contrary, one 
can keep the coupling fixed and increase the system tem- 
perature, "switching on" more and more bath oscillators 
until their number reaches the maximum at the Debye 
temperature To = hiOrnax/kB- In such a case not only 
the T-mode width increases due to higher dissipation, 
but also the peak position shifts with temperature. 

The temperature dependence of the T-mode frequency 
and its FWHM are shown in Figs. [2]|4l In the sub- and 
super-Ohmic regimes we used the power indexes n — 0.03 
and n = 1.97 (which slightly deviate from the corre- 
sponding integer values and 2) to eliminate the lower 
bound cut-off frequency ujq. The others fitting para- 
meters (the vibrational frequency J7 and the coupling 
strength rjc) are chosen to reproduce the experimental 
data g. 



It is seen from Fig.[2]that in a strong sub-Ohmic regime 
the T-mode frequency is Wwe-shifted with T up to tem- 
peratures of about hfl/ks, where the T-mode frequency 
becomes a descending function of temperature. 



0)r(T)/m, 



A(T)/m 



(4.8) 




FIG. 2: Left panel: temperature dependence of the T-mode 
frequency in the sub-Ohmic regime with n = 0.03, dimen- 
sionless vibrational frequency Q,/ujmax = 0.2 and coupling 
constant i]c ~ 0.002. Right panel: temperature dependence 
of the T-mode FWHM at the same parameters. 

The T-mode FWHM is found to be a nonlinear func- 
tion of T at all temperatures studied. At ksT /tvjjmax ~ 
0.4 the T-mode becomes very broad in spite of very small 
coupling, and almost disappears from the spectrum func- 
tion Re[C'i,„((u;)] (see next Section for analytic estima- 
tions). 

To explain the T-mode Wue-shift at low-to-moderate 
temperatures let us recall that in the sub-Ohmic regime 
with n = there are residual correlations of the fiicker 
noise [27'| . The correlation time T^or of the thermal bath 
random forces is much larger than the relaxation time Ty 
of the velocity autocorrelation functions (see Eqs. (14.41) . 
(|4.8p '). Thus, a non-relaxed energy of the lattice is be- 
ing delivered to the adparticle increasing its effective vi- 
brational energy until the inverse process begins at the 
temperatures of about HQ/ks- 



01^(1)/, 




A(T)/, 




FIG. 3: Left panel: temperature dependence of the T-mode 
frequency in the super-Ohmic regime with n — 1.97, dimen- 
sionless vibrational frequency ^/ujmax = 0.2 and coupling 
constant i]c ~ 0.166. Right panel: temperature dependence 
of the T-mode FWHM at the same parameters. 

In Fig. [3] we present the results for the T-mode tempe- 
rature induced shift and broadening in the strong super- 
Ohmic regime. There are ultra-short range time corre- 
lations Tcor of the thermal bath random forces at super- 
Ohmic coupling j23| . and much longer times t^ at which 
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the velocity of the adparticle decays (see Eqs. (|4.5I) . 
(|4.8p ). In such a case the excess energy of the adsorbate is 
being transferred to the lattice, yielding a decrease of the 
effective vibrational energy of the adparticle with tem- 
perature and, consequently, the red-shift of the T-mode 
frequency. The FWHM of the T-mode accommodates 
to almost a linear function of temperature after a short 
transition regime at low T . 

At the Ohmic coupling (see Fig. HJ the T-mode is 
slightly Wwe-shifted with temperature, and its FWHM 
is a linear function of T like in the super-Ohmic case. At 
the Ohmic coupling the difference between two typical 
times Tcor and t„ is smaller than in the other regimes. 
As a result, the T-mode shift is also smaller than in the 
sub- or super-Ohmic regimes (at the same values of the 
parameters of the model), which will be shown in the 
next Section. Such a behavior is consistent with the ex- 
perimental data of Ref. [29] where the e-h pair contribu- 
tion (with a weak Ohmic coupling [2l|) was found to be 
almost independent of temperature. 



quantum Zeno and anti-Zeno effects, however, without 
taking into account the tunnelling-mediated surface dif- 
fusion and memory effects. 

In our paper we consider a more general case when the 
"driving force" is dealt with slow underbarrier sitc-to-site 
hopping of the adsorbate and should be treated on equal 
footing with other interaction mechanisms. The interac- 
tion term (|2.3p of Hamiltonian is not of a bilinear form 
in creation/annihilation operators. It takes into account 
coupling both to density and oscillation modes of the ad- 
sorbate that corresponds to more realistic systems [2, 18j. 
Neither a power series expansion for the spectral weight 
functions at vibration frequency $7 nor any approxima- 
tion for the kinetic kernels (like the Wigner-WeisskopfF 
one) are made in our paper. As a result, we manage to 
describe the T-mode behavior quite accurately using only 
one fitting parameter - the dimensionless coupling con- 
stant. It is shown in the next Section that its value can 
be extracted from the residual [5] shift and broadening 
of the T-mode at zero temperature. 



0.14161 '^(T)'". 




A(T)/m 
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FIG. 4: Left panel: temperature dependence of the T-mode 
frequency in the Ohmic regime with dimensionless vibra- 
tional frequency Q,/uJmax = 0.141 and coupling constant 
77c = 0.0066. Right panel: temperature dependence of the 
T-mode FWHM at the same parameters. 



V. AN ANALYTICAL EVALUATION OF THE 
T-MODE TEMPERATURE DEPENDENCE 

Analytical results can be obtained in the zero tempera- 
ture and high temperature limits. In the first case, a ther- 
mal population factor from Eq. (|2.1ip is eliminated. On 
the other hand, at high temperatures the thermal popu- 
lation factor can be approximated as coth.{tujj /2kBT) ~ 
2kBT/hu!. Then using the long time asymptotics (|4.4p - 
(|4.6p for the kinetic kernels instead of the exact expres- 
sions (|2.9|) we can obtain analytical expressions for the 
T-mode frequency and FWHM. 



Summing up this Section, we would like to discuss 
the difference between our approach and the one used 
in Ref. [5| when studying the temperature behavior of 
the adsorbate T-mode. In Ref. [5] two fitting parameters 
were used, which characterize the behavior of the spectral 
weight function Ja;(w) times a thermal population factor 
in the vicinity of the T-mode frequency. These param- 
eters were found to depend very weakly on the surface 
temperature for the three systems studied. Moreover, the 
dynamic structure factor itself has been calculated in the 
Wigner-WeisskopfF approximation [ij] . That introduced 
an additional inaccuracy in the description of the sys- 
tem dynamics at intermediate times, when the memory 
effects have to be taken into account rigorously [19, 20]. 

Even though in Ref. |5| the quantum model for the ad- 
sorbate motion has been proposed, further calculations 
were carried out in the impulsive collision approximation 
without reference to the driving Hamiltonian term. The 
influence of non-selective measurements (e.g. collisions) 
on the evolution of quantum system was considered more 
precisely in [30j . where the short- and long-time dynam- 
ics of the adsorbate has been studied in the context of 



A. Zero temperature limit 

In the Ohmic regime and in the weak coupling limit 
77c ^ 1 the following expression for the T-mode frequency 
is valid 



CJTiO) = V'fi2 - A(0)2/4, 



(5.1) 



that coincides with the expression for the frequency of the 
damped harmonic oscillator. The expression (|5.ip for the 
T-mode localization frequency defines a shift that can be 
referred to as the residual one, insofar it is present even 
at K. The corresponding residual (or intrinsic) FWHM 
in the sub- and Ohmic regimes (ri < 1) is scaled as 



A(0) ^ 77cr!" 



(5.2) 



The damping coefficient (|5.2p is not zero in the T -> 
limit, because the excited adsorbate can transfer its en- 
ergy to the lattice, induce excitations from the zero-point 
motions of the lattice modes, or create e-h excitations in 
the electron distribution. 



Let us note that the scahng law (I5.2p for the T-mode 
FWHM agrees completely with the results of Ref. [5J] ob- 
tained in the zero temperature limit. Besides, a practical 
significance of such results becomes quite clear: inspect- 
ing the T-mode localization, one can draw a conclusion 
about the vibrational frequency J7 whereas FWHM in- 
cludes the information about the "adsorbate-substrate" 
coupling strength rjc- 

Our model does not admit a straight zero tempera- 
ture limit in the super-Ohmic regime because there is no 
polaron band narrowing in such a case [iS], and a new 
derivation of the coherent term of the generalized diffu- 
sion coefficient becomes indispensable. Nevertheless, the 
obtained results are quite reliable to be used for the anal- 
ysis of the low temperature T-mode features from both 
theoretical and experimental viewpoints. 



B. High temperature limit 

1. Sub-Ohmic regime with n — 

Taking into account the Gaussian form (|4.4p of the long 
time asymptotics for the kinetic kernels and Eq. (|4.7p for 
the frequency dependence of the T-mode it is straightfor- 
ward to obtain the expansion for the T-mode frequency 

/^x n I ^/cr|lna;o| /^ 5 r/cTl lna;o| \ 
u:r{T) = n+ |1-^ ^, I (5.3) 

as a series in the coupling strength rjc- The corresponding 
FWHM A(T) can be written down as 



A(T) = 



2^T]cT\\nujo\ 



(5.4) 



it does not admit the r]c series expansion. The T-mode 
FWHM is found to be a nonlinear function of tempera- 
ture, which agrees completely with the results presented 
in Fig. El 

Analyzing the T-mode temperature dependence (|5.3p 
it is possible to find the temperature 



rp^ 



02 



Sr^dlnwol 
at which a maximal shift of the T-mode 

Awt(T*) = ^t{T*) - O = ^O 



(5.5) 



(5.6) 



takes place. The shift (j5.6p is quite large, and the T-mode 
FWHM at this temperature 

A(r*) = y|cxp(^-^)o«io (5.7) 

approaches the vibrational frequency, leading to a rapid 
decay of the excitation as it is clearly seen in Fig. [21 



2. Ohmic regime 

Analysis of the T-mode temperature dependence at the 
Ohmic coupling shows that at reasonable values of the 
vibrational frequency O the T-mode shift 



Aujt{T) 



20 



(5.8) 



is smaller than in sub- or super-Ohmic systems, since 
the expansion in (j5.8p starts from a quadratic term in 
the coupling strength. The FWHM is governed by the 
expression 



A(T) - 77cTO"- 



(5.9) 



which is also valid for the super-Ohmic systems with f < 
n < 2. 

It should be stressed that the temperature dependence 
of the T-mode frequency (|5.8p differs from the results 
presented in Fig. S] The T-mode is red-shifted in accor- 
dance with Eq. (|5.8p . whereas the ft/we-shift is observed 
in Fig. 21 An explanation of this contradiction is quite 
simple: we used the exponential form (|4.6p of the ki- 
netic kernel when calculating (j5.8p . while at the numeri- 
cal evaluation we used the exact expression (j2.9p . which 
up to 5-10% differs from asymptotic value (|4.6I) at short 
times. This deviation can be considered negligible when 
calculating the temperature dependence of the diffusion 
coefficients [1^, [201 ■ But it is not true when one evaluates 
the temperature dependence of one-particle excitation of 
the adsorbate, that once more shows the exceptional sen- 
sitivity of the T-mode behavior to any kind of approxi- 
mations. 



3. Super-Ohmic regime with n = 2 

At the super-Ohmic regime with n = 2, weak coup- 
ling limit rjc <S^ 1, and moderate vibrational frequencies 
O < LOmaxi the temperature dependence of the T-mode 
localization can be presented as 



ujt{T) = + r7cr01n(20). 



(5.10) 



It is seen from the last expression that the T-mode 
is red-shifted with temperature up to the vibrational 
frequencies ^/tOmax — 1/2, being at the same time a 
strongly nonlinear function of O. Exact numerical eval- 
uations show that the a;T(T)-curve at first reaches a 
plateau as the vibrational frequency increases and then 
becomes a decreasing function of temperature at O ^ 

^maxj ^• 

On the other hand, the T-mode FWHM is very small 
even at the Debye temperature 



A{kBT/hw^ 



1) -- rjcVL < O. 



(5.11) 



Thus, the T-mode is not broadened too much and should 
be well defined. It is also to be pointed out that 
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in the super-Ohmic regime the T-mode shift is much 
smaller than that in the sub-Ohmic case, since Jlj lnJ7| ^ 
llnwol/ri. 

Summarizing, the analytic investigation of the T-mode 
behavior confirms completely the main tendencies ob- 
served in Figs. [2EI Namely: i) a linear increase of 
the T-mode FWHM with temperature in the Ohmic and 
super-Ohmic regimes and a non-linear increase of A(T) 
in the sub-Ohmic regime; ii) the corresponding red- or 
blue-shiits of the T-mode location; iii) a considerable T- 
mode broadening in the sub-Ohmic regime even at weak 
"adsorbate-substrate" coupling and a small broadening 
at 1 < n < 2. 

Since there is a one-to-one correspondence between the 
system ohmicity and the kind of the coupling, one can 
deduce about a nature of the "adsorbate-substrate" in- 
teraction from the temperature induced shift of the T- 
mode. When several kinds of the "adsorbate-substrate" 
interaction act simultaneously, the data obtained could 
be more complicated than those reported in this Section. 
However, such a situation is rather infrequent since usu- 
ally the typical relaxation times of different interactions 
are separated by several orders of magnitude @, [a] , and 
it is possible to talk about a single effective interaction. 



VI. T-MODE AND MULTIPLE ADSORBATE 
HOPPING 



Investigation of the adparticle diffusion on the basis 
of the Klein-Kramers equation [l|, [l^, [ij] allows us to 
look at the problem of T-mode from another viewpoint, 
considering this one-particle excitation as a precursor of 
multiple or long jumps of the adsorbate. Such a jumping 
regime should not be confused with the situation, when at 
very low barriers the particle performs a quasi-continuous 
hopping at several sites [l|, [ij] , and a diffusive stage of 
evolution is not reached. One rather has to talk about 
the situation when the velocity of the particle, which has 
localized at a certain adsorption site s, is not thermal- 
ized yet, and the adsorbate jumps to the next nearest 
neighboring site s' very soon. Thus, we will use the term 
"multiple" rather than "long" hopping of the adsorbate. 

In Ref. [141 two conditions for the multiple jumps of the 
adsorbate at the thermal activated diffusion have been 
established, formulated as strong inequalities 

Tth < T„, Tosc < Ty (6.1) 

for three typical times of the system dynamics 



1 

— , Tth = a 

Vc 



knT' 



(6.2) 



where Tqsc rneans the period of oscillation of the adpar- 
ticle with mass m at the bottom of the potential well of 
the depth U; Ty is a relaxation time of the velocity auto- 
correlation function, and Tth denotes the time taken by 
the particle to cross over a lattice spacing a with a mean 
thermal velocity vth- 



It is seen from the presented above analysis that the 
T-mode is defined by the second inequality of (|6.ip . If 
the system parameters allow the typical times to obey 
the first inequality too, then the multiple jumps can be 
observed in the system. 

Let us suppose that conditions (|6.ip are fulfilled. Then 
during both the jump to the nearest neighboring site (left 
inequality) and oscillation within a certain adsorption 
site (right inequality) the particle does not have enough 
time to be thermalized, and the multiple hopping sce- 
nario of the diffusion can be realized. In such a case, 
the Einstein-Smoluchowski equation is found to be in- 
sufficient pj, |17]] for the description of the adparticle dif- 
fusion, and one has to use the Klein-Kramers equation 
for the distribution function, which depends both on the 
velocity and on the coordinate of the adsorbate. 

In the case of quantum diffusion the multiple hopping 
conditions (j6.ll) remain valid too, but the values Tqsc and 
Tth from Eq. (|6.2I) should be substituted by the typical 
times 



2n 



2nh 



Ttv 



(6.3) 



where Tqsc denotes the inverse vibrational frequency, and 
Ttun is a time taken by the particle to perform an under- 
barrier hopping to the nearest neighboring site. 

It should be pointed out that Refs. [l|, [I4I were di- 
rected mainly to the investigation of inelastic peak of the 
dynamic structure factor (or, alternatively, the side peak 
of velocity autocorrelation function) at various coupling 
constants and barrier heights, and no study of the tem- 
perature dependence of the inelastic peak was performed 
to attribute it to the T-mode features. 

On the other hand, a generalization of the Klein- 
Kramers equation to description of the quantum diffu- 
sion (especially in the weak coupling limit, when there 
is the energy-diffusion-controlled regime from the view- 
point of transition state theory [l5|) is an intricate prob- 
lem, which, up to our knowledge, has not been solved yet. 
In Ref. [20] we followed an alternative way, having stud- 
ied a behavior of the "velocity-velocity" time correlation 
functions, defined on the adjacent sites, with taking into 
account the memory effects. It has been shown that at 
given parameters of the model the multiple jumps are 
absent (the relation Ttun '^ Ty was found to be valid), 
but they should appear when the tunnelling amplitude 
increases or the coupling strength decreases. The most 
interesting question - about the influence of temperature 
on the multiple hopping of adsorbate - can be answered 
only after an additional analysis of the T-mode behavior. 



VII. CONCLUDING REMARKS 

In this paper, we carry out numerical and analytical in- 
vestigation of the temperature dependence of the T-mode 
of adsorbate which diffuses from one adsorption site to 
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another by a tunnelling mechanism. The T-mode charac- 
teristics (its localization frequenc y an d width) are known 
to depend on temperature 0, H [20| • Having adopted a 
concept that the T-mode appearance is dealt with mem- 
ory effects of the adsorbate motion, we propose an expla- 
nation of cither {red or blue) temperature induced shifts 
of the above mentioned one-particle excitation on the ba- 
sis of comparative analysis of two typical timescales: the 
correlation time Tcor of the thermal bath random forces 
and the relaxation time t„ of the velocity autocorrelation 
functions. 

We show that in the strong sub-Ohmic regime with 
n <C 1 the relation r^ ^ Tcor is valid, leading to a blue- 
shift of the T-mode with temperature, since the excess 
energy of the lattice is being delivered to the adparticle, 
thereby increasing its effective vibrational energy. On the 
contrary, in the strong super-Ohmic regime with n w 2 
the inequality r^ 3> Tcor becomes valid. At the timescales 
of about Tcor the adparticle does not have enough time 
to be thermalized and starts to transfer its energy to the 
lattice, yielding the observed red-shift of the T-mode. 

We have also analyzed in detail the temperature behav- 
ior of the T-mode width, which changes from a strongly 
nonlinear function of T (in the case of strong sub-Ohmic 
regime with n <C 1) to the linear dependence in temper- 
ature (for Ohmic and super-Ohmic systems). 

The ohmicity type is closely related to the nature of the 
"adsorbate-substrate" interaction: the case with n = 
corresponds to the flicker noise, the Ohmic regime is be- 
ing realized in the systems with "electronic" friction or 
at the coupling with nonlinear phonons, while the super- 
Ohmic case with n = 2 corresponds to the linear De- 
bye phonon coupling. Therefore, one can draw a conclu- 
sion about the interaction mechanism from analysis of 
the temperature induced shift of the T-mode. On the 
other hand, extrapolating the obtained data to the zero 
temperature limit, it is possible to evaluate the micro- 
scopic parameters (vibrational frequency fl and coupling 
strength rjc), which govern behavior of the "adsorbate- 
substrate" system. 

We consider a slow underbarrier site-to-sitc hopping 
of the adparticle to be a "driving force" of the surface 
motion of the adsorbate. Unlike the impulsive collision 
approximation [5| , a "driving term" (12. 2p of the Hamil- 
tonian is treated on equal footing with other interaction 



mechanisms and defines a time evolution of the system. 
Since there is a lack of data about the T-mode behavior 
for the case, when the diffusion is mediated by tunnelling 
and the non-Markovian effects play an important role [3Q] 
(or, what is almost the same, when different regimes of 
the ohmicity manifest themselves), we believe that our 
investigations have a good perspective from a viewpoint 
of the study of transition regimes of the adsorbate. 

It has to be emphasized that a similar approach should 
be quite promising in the case of thermally activated dif- 
fusion as well. It was shown in Ref. [3l| that in non- 
Ohmic systems with < n < 2, n ^ 1, the thermal 
bath random forces were correlated as {Fr{t)Fr{0)) ~ 
^'"^^ , whereas the long time asymptotics of the veloc- 
ity autocorrelation functions was scaled as {v{t)v{0)) ^ 
^"'^ -p^TT- Thus, the relations r„ <^ Tcor or r^, ^ Tcor 
remain valid also at thermally activated diffusion, allow- 
ing one to perform a similar comparative analysis of the 
typical timescales and to study the temperature behavior 
of the T-mode. 

Finally, we relate the T-mode to the onset of the mul- 
tihops of the adsorbate. If an eventual temperature in- 
duced shift of the T-mode toward lower frequencies and 
its broadening are not too large (quasi- and inelastic 
peaks of the dynamic structure factor are well resolved) , 
one can state that there is a portion of multihops in 
the system. The computer simulations within the Monte 
Carlo wave function formalism [32] and direct evaluation 
of the "flux-flux" time correlation functions [3^, [SJ] sup- 
port this suggestion. Thus, investigation of the T-mode 
behavior along with a study of the "velocity-velocity" 
cross correlation functions can shed more light upon the 
nature of such an interesting phenomenon of the adsor- 
bate motion as a diffusion by the multi-hopping scenario. 
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